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We calculate the quarkonium dissociation rate in the P and CP-odd domains (bubbles) 
that were possibly created in heavy-ion collisions. In the presence of the magnetic field 
produced by the valence quarks of colliding ions, parity-odd domains generate electric field. 
Quarkonium dissociation is the result of quantum tunneling of quark or antiquark through 
the potential barrier in this electric field. The strength of the electric field in the quarkonium 
comoving frame depends on the quarkonium velocity with respect to the background mag- 
netic field. We investigate momentum, electric field strength and azimuthal dependence of 
the dissociation rate. Azimuthal distribution of quarkonia surviving in the electromagnetic 
field is strongly anisotropic; the form of anisotropy depends on the relation between the 
electric and magnetic fields and quarkonium momentum P±. These features can be used to 
explore the properties of the electromagnetic field created in heavy ion collisions. 



Solid theoretical arguments pQ and numerical calculations [3] indicate a possible existence 
of very strong magnetic fields in heavy-ion collisions. Electromagnetic fields of such enormous 
intensity have never been experimentally studied and therefore present a great interest, which 
extends far beyond applications in the nuclear physics. What are the possible manifestations of 
such magnetic field? An effect that has recently attracted a lot of attention is the Chiral Magnetic 
Effect (CME) [H |1HZ!- If a metastable P and CP-odd bubble is induced by axial anomaly in the 
hot nuclear matter, then in the presence of external magnetic field Bq the bubble generates an 
electric field which is parallel to the magnetic one. According to [5] the value of the electric field 
Eq in the bubble is 



where the sum runs over quark flavors / and we assumed that only three lightest flavors contribute. 
The value of the 0-angle fluctuates from event to event. CME refers to the macroscopic manifes- 
tation of this effect - separation of electric charges with respect to the reaction plane. This effect 
is a possible explanation of experimentally observed charge asymmetry fluctuations [8l-ll0|. 
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INTRODUCTION 
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Other effects of the magnetic field that do not require existence of the parity-odd bubbles and 
have been recently discussed are: synchrotron radiation by fast fermions [TT], polarization of the 
fermion spectra (also known as the Sokolov-Ternov effect [12])) enhanced dilepton production 
|13j and azimuthal anisotropy of the quark-gluon plasma (QGP) |14| 115]. These possible effects 
await their experimental study. Thus, strong magnetic field must have a powerful impact on the 
behavior of the quark-gluon plasma. 

Since CME cannot be the only effect induced by the magnetic field we are motivated to look 
for magnetic field manifestations in other observables. We have recently pointed out in Ref. [16] 
that the Lorentz ionization, i.e. dissociation of a moving bound state in external magnetic field, 
is phenomenologically significant in relativistic heavy-ion collisions in the presence of strong [U 
[2] quasi-static [11] magnetic field generated by colliding ions. The ionization or, equivalently, 
dissociation happens because quarkonium constituents have finite probability to tunnel through 
the potential barrier in the presence of electric field, which appears upon boosting to the comoving 
frame. The dissociation rate depends on the magnetic field strength, bound state velocity and its 
binding energy. The dissociation rate indirectly depends on the properties of the nuclear matter 
by the way of dependence of the binding energy on temperature. Since at higher temperatures the 
binding energy is smaller, the dissociation rate is higher. Still, the Lorentz ionization can happen 
even if no matter is formed, provided that the magnetic field is strong enough and/or the bound 
state is fast enough. 

In this paper we address a different type of ionization which is possible only if a metastable P 
and CP-odd bubbles are formed and trigger emergence of electric field Eq given by Due to 
this field, there is a finite dissociation rate of quarkonium at rest in the laboratory frame. When 
boosted to the comoving frame quarkonium dissociation is a combined effect of the boosted electric 
Eq and magnetic Bq fields. In the comoving frame, electric and magnetic fields are directed at 
some angle with respect to each other depending on the quarkonium kinematics. The main goal of 
this paper is to investigate the dissociation rate in this case. 

II. DISSOCIATION RATE 

Ionization probability of quarkonium equals its tunneling probability through the potential 
barrier. In the WKB approximation the later is given by the transmission coefficient and was 
calculated in |16| . In this method contribution of the quark spin can be easily taken into account. 
Another method of calculating the ionization probability, the imaginary time method \17\ 118] . 
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was employed in |19H21| by Popov, Karnakov and Mur. In particular, they derived in the non- 
relativistic approximation the pre-exponential factor that appears due to the deviation of the quark 
wave function from the quasi-classical approximation. Such a calculation requires matching quark 
wave function inside and outside the potential barrier [23J. Extension of this approach to the 
relativistic case is challenging due to analytical difficulties of the relativistic two-body problem. 
Fortunately, it was argued in TT\ that the non-relativistic approximation provides a very good 
accuracy in the e^^m region which is relevant in the quarkiononium dissociation problem. 

Magnetic field produced in heavy-ion collisions may have a complicated spatial and temporal 
structure (see e.g. [23 )• However, we will assume that the field is constant on the scales relevant for 
the problem of J/ip dissociation. This approximation is supported by our recent arguments [II] that 
the relaxation time of magnetic field is of the order of a few fm/c due to high electrical conductivity 
of the QGP*. Also, spatial inhomogeneity in the transverse plane reveals itself at distances of the 
order of the nucleon size and perhaps even larger in view of uniformity of the matter distribution 
in the nuclei. A more quantitative estimate of the role of spatial and temporal dependence of the 
magnetic field on J/ip dissociation requires numerical solution of magneto- hydrodynamic equations 
and is beyond the scope of this paper. 

Given the electromagnetic field in the laboratory frame -Bo, Eq, the electromagnetic field B, E 
in the comoving frame moving with velocity V is given by 

E =E irUbo + p^V x 6q) - ( 7 L - 1)^^^} (2a) 

B =B | 7L (&o - PoV x b ) - (lL ~ 1)^^^} (2b) 

where &o = Eq/Bq is a unit vector in the magnetic field direction, po = Eq/Bq = 2a|G|/37r (see 



lk) and 7l = 1/Vl - V 2 . It follows from <m that 



E = E Jl + 72 (b x V) 2 (l + po 2 ) (3a) 



B = B ^l + 1 l{b Q xVf{l + pl) (3b) 
Using 43) we find that the angle 6 between the electric and magnetic field in the comoving frame 



is 



' -' '[1 + ll {bo x V) 2 (l + p„- 2 )][l + -ri (to x Vf(\ + pl)\ 



* Calculations in Refs. [IH3] yield very short relaxation time of magnetic field because they neglect the electromag- 
netic response of the quark-gluon medium, see [11] for details. 
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where we used the relativistic invariance of E ■ B. 

It is useful to introduce dimensionless parameters 7, e and p as [21 J 

1 /2ib E eE ( m\ 3/2 

p \ m B m 2 \ 2sb J 

where m is quark mass and £& is quarkonium binding energy. We will treat the quarkonium 
binding potential in the non-relativistic approximation, which provides a very good accuracy to 
the dissociation rate 121] . The quarkonium dissociation rate in the comoving frame in the 
non-relativistic approximation is given by |19j 

w = ^ j p( 7 ^)C 2 (7,0)e-^ 9(7 ' e) (6) 



where function g reads 



3t 
27 



1 /r 2 \ 1/2 r 2 
-f^-1] sinfl-^cos 2 . 
7 V7 / %T 



(7) 



and functions P and C are given in the Appendix. The contribution of quark spin is taken into 
account by replacing e& — > e' b = — ■ B [16]. Function g represents the leading quasi- classical 
exponent, P is the pre- factor for the S- wave state of quarkonium and C accounts for the Coulomb 
interaction between the valence quarks. Parameter tq satisfies the following equation 

r 2 - sin 2 (9(r coth r - l) 2 = 7 2 (8) 

which establishes its dependence on 9 and 7. Note, that in the limit E — > the dissociation rate 
Q exponentially vanishes. This is because pure magnetic field cannot force a charge to tunnel 
through a potential barrier. 

Eq. ^ gives the quarkonium dissociation rate in a bubble with a given value of 0. Its derivation 
assumes that the dissociation process happens entirely inside a bubble and that is constant inside 
the bubble. Since in a relativistic heavy ion collision many bubbles can be produced with a certain 
distribution of 0's (with average (0) = 0) more than one bubble can affect the dissociation process. 
This will result in a distractive interference leading to reduction of the CP-odd effect on quarkonium 
dissociation. However, if a typical bubble size Ro is much larger than the size of quarkonium Rj, 
then the dissociation is affected by one bubble at a time independently of others, and hence the 
interference effect can be neglected. In this case Q provides, upon a proper average, a reasonable 
estimate of quarkonium dissociation in a heavy ion collision. We can estimate the bubble size as 
the size of the sphaleron, which is of the order of the chromo-magnetic screening length ~ l/g 2 T, 
whereas the quarkonium size is of the order a s /eb- Consequently, at small coupling and below the 



zero-field dissociation temperature (i.e. when e& is not too small) Rq is parametrically much larger 
than Rj. A more quantitative estimate of the sphaleron size is Rq ~ 1.2/ a s N c T ~ 0.4 fm [22] : 
whereas for J/tp Rj ~ a s /eb ~ 0.1 — 0.2 fm. Thus, based on this estimate bubble interference can 
be neglected in the first approximation. However, since the ratio Rj j 1 Rq is actually not so small 
this effect nevertheless warrants further investigation. 

To obtain the experimentally observed J/tp dissociation rate we need to average ^ over the 
bubbles produced in a given event and then over all events. To this end it is important to note 
that because the dissociation rate depends only on p^ it is insensitive to the sign of the E field or, 
in other words, it depends only on absolute value of O but not on its sign. Therefore, it stands to 
reason that although the precise distribution of G's is not known, ([6]) gives an approximate event 
average with parameter representing a characteristic absolute value of the theta-angle. 



III. LIMITING CASES 



Before we proceed with the numerical calculations, let us consider for illustration several limiting 
cases. If quarkonium moves with non-relativistic velocity, then in the comoving frame electric and 
magnetic fields are approximately parallel 8 ~ 0, whereas in the ultra-relativistic case they are 
orthogonal 6 ~ vr/2, see Q. In the later case the electromagnetic field in the comoving frame 
does not depend on Eq as seen in ^ and therefore the dissociation rate becomes insensitive to 
the CME. In our estimates we will assume that pq < 1 which is the relevant phenomenological 
situation. Indeed, it was proposed in [5] that po ~ a <C 1 produces charge fluctuations with 
respect to the reaction plane of the magnitude consistent with experimental data. 

1) ^ ^ 0) i- e - electric and magnetic fields are approximately parallel. This situation is realized 
in the following two cases, (i) Non-relativistic quarkonium velocities: V -C po or (ii) motion of 
quarkonium at small angle (p to the direction of the magnetic field &o : 4> *C Po/tl^- In both cases 
E ~ Eq and B ~ Bq. This is precisely the case where the dissociation rate exhibits its strongest 
sensitivity to the strength of the electric field Eq generated by the local parity violating QCD 
effects. Depending on the value of the 7 parameter defined in ^ we can distinguish the case of 
strong electric field 7 > 1 and weak electric field 7 C 1 |20j . In the former case, g = (3/8)7, 
P = (8/e) 1 / 2 7e _ ' y2 / 2 , C = e 71 " 7 / 2 /j. Substituting into Q the dissociation rate reads 
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In the later case, g = P = C = 1 and 



8e b 2 8e b m 2 f2e b \ 3/2 -^(^) 3/2 



e ei^o \ m J 

where the electromagnetic field in the comoving frame equals one in the laboratory frame as was 
mentioned before. 

2) 6 ~ vr/2, i.e. electric and magnetic fields are approximately orthogonal. ^ This occurs for an 
ultra-relativistic motion of quarkonium V — > 1. In this case 



B = E = B o7L \b xV\^l + p 2 (11) 

This case was discussed in detail in our previous paper [IB] . In particular for 7 <C 1 we get 

8e t m 2 (2e t \ 3 ' 2 -isil'^V" 



—- s-U?J e (12) 



Due to ([9]) and (12) dependence of if on Eq is weak unless po S> 1. 



IV. J/V> DISSOCIATION RATE 

One of the most interesting applications of the formalism described in the previous sections 
is calculation of the dissociation rate of J/ ifi which is considered a litmus test of the quark-gluon 
plasma f25j. 

Let z be the heavy ions collision axis; heavy-ion collision geometry implies that bo ■ z = 0. The 
plane containing z-axis and perpendicular to the magnetic field direction is the reaction plane. We 
have 

(60 x V) 2 = V% + Vl sin 2 <f> (13) 

where <p is the angle between the directions of Bq and VI and we denoted vector components in 
the xy-plane by the subscript JL We can express the components of the quarkonium velocity V in 
terms of the rapidity rj as V z = tanhr;, V± = P±/(M± cosh?]), where P and M are the quarkonium 
momentum and mass and M"j_ = M 2 + P^_ . 

Results of numerical calculations are exhibited in Figs. [l}j3} In Fig. [T] we show the dissociation 
rate of J/ip for several values of the electric field Eq induced by the Chiral Magnetic Effect. 
Note, that the typical size of the medium traversed by a quarkonium is a few fm. Therefore, 

t Note, that the limit 7 3> 1 is different in 9 = 5 and 6 < f cases |19] , 
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FIG. 1: Dissociation rate of J/tp at eBo = 15m£, <fi — ir/2 (in the reaction plane), rj = (midrapidity) as a 
function of (a) P± at = 0.16 GeV and (b) at P± = 1 GeV. 

w ~ 0.3 — 0.5 fm _1 corresponds to complete destruction of J/^'s. This means that in the magnetic 
field of strength cBq ~ 15m^ all J/0's with Pj_ > 0.5 GeV are destroyed independently of the 
strength of Eq. This P± is lower than we estimated previously in |16| neglecting the pre- factors 
in the dissociation rate. Since magnetic field strength decreases towards the QGP periphery, most 
of J/ip surviving at later times originate from that region. Effect of the electric field Eq of the 
parity-odd bubble is strongest at low P±, which is consistent with our discussion in the previous 
section. The dissociation rate at low P± exponentially decreases with decrease of -E^. Probability 
of quarkonium ionization by the fields below Eq < O.IBq (i.e. pq < 0.1) is exponentially small. 
This is an order of magnitude higher than the estimate pq ~ a proposed in io] . 

As the plasma temperature varies, so is the binding energy of quarkonium although the precise 
form of the function £b{T) is model-dependent. The dissociation rate picks at some e® < (see 
Fig. [l|b)), where e™ c is the binding energy in vacuum, indicating that J/ifi breaks down even 
before e\, drops to zero, which is the case at Bq = 0. This is a strong function of Eq as can be 
seen in Fig. [2j It satisfies the equation dw/d£b = 0. In the case 7 1 (10) and (12) imply that 

2/3 



-0 



m ( 5eE\' 



7 < 1 (14) 
At 7 » 1 and 6 = tt/2 we employ S to derive the condition (e°) 2 + ei?e°/2m — eE 2 /B = 0. In 



view of (11) E ~ B and we obtain 




+ 1-1 



'eB, 7 > 1 



(15) 



where in the last step we used that eB <C m 2 . For a given function £b(T) one can convert e® into 
the dissociation temperature, which is an important phenomenological parameter. 
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FIG. 2: Contour plot of the dissociation rate of J/ip as a function of e& and eEg at eBo — 15m„, <j> = tt/2 
(in the reaction plane), rj = (midrapidity) and P± =0.1 GcV. Numbers inside boxes indicate the values 
of w in 1/fm. 
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FIG. 3: (a) Angular distribution of J/ip dissociation rate at eBo — 15m^, i] = at different Eg and P± 
(in GeV's). Magnetic field Bo points in the positive vertical direction. Reaction plane coincides with the 
horizontal plane, (b) Rescaled second Fourier-harmonic vi of the azimuthal distribution as a function of 
P± . (P) is the azimuthal average of the survival probability and t is the time spent by J/tp in the P-odd 
bubble. 



In the absence of the CME the dissociation probability peaks in the direction perpendicular 
to the direction of magnetic field bo, i.e. in the reaction plane. Dissociation rate vanishes in the 
&o direction. Indeed, for V • &o = Q implies that E = 0. This feature is seen in the left 
panel of Fig. [3| At finite Eq the dissociation probability is finite in the bo direction making the 
azimuthal distribution more symmetric. The shape of the azimuthal distribution strongly depends 
on quarkonium velocity: while at low V the strongest dissociation is in the direction of the reaction 
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plane, at higher V the maximum shifts towards small angles around the bo direction. Extrema of 
the azimuthal distribution are roots of the equation dw/d(p = 0. At 7 <C 1 it yields minimum at 
4>o = 0, maximum at (f)Q = ir/2 and another maximum that satisfies the condition (neglecting the 
spin-dependence of £b) 



2m 2 



2 /9c, \ 3/2 



eJ EoVl + 7i(K 2 + ^sm 2 0o)(l + Po") = ^" [^) (16) 



In order to satisfy ( |16| 4>o must decrease when V increases and visa versa. This features are seen 
in the left panel of Fig. |3| 

Spectrum of quarkonia surviving in the electromagnetic field is proportional to the survival 
probability P = 1 — wt, where t is the time spent by the quarkonium in the field. Consider P as 
a function of the angle \ between the quarkonium velocity and the reaction plane \ = k/2 — <f>. 
Fourier expansion of P in \ reads 

P(X) = l P o + jt Pn cos ( n *) ' p n = ~ [ P(X) cos(nx) d X . (17) 

n=l n 

Ellipticity of the distribution is characterized by the "elliptic flow" coefficient vi defined as 

V2= Wo = — «(p) — = ~W)L W cos2xdx (18) 

where (P) denotes average of P over the azimuthal angle. These formulas are applicable only as 
long as wt < 1 because otherwise there are no surviving quarkonia. In the right panel of Fig. [3] we 
show V2 (P) /t, which is independent of t, as a function of P±. As expected, in the absence of the 
CME, V2 is negative at low P± and positive at high P±. V2 changes sign at P± that depends on 
the strength of the electric field. It decreases as Eq increases until at Eq ~ Bq it becomes positive 
at all P±. Fig. (3^b) provides the low bound for V2 because (P) < 1 and t > 1 fm . We thus expect 
that magnetic field strongly modifies the azimuthal distribution of the produced J/^'s. Role of the 
magnetic field in generation of azimuthal anisotropics in heavy-ion collisions has been pointed out 
before in 



V. SUMMARY 



In this paper we studied the effect of the parallel electric and magnetic fields on the dissociation 
rate of quarkonia, and particularly of J/ip. Our main observation is that the CME effect on the 
dissociation rate is significantly different than the effect of the pure magnetic field if Eq > 0.1 -Bo 
that implies an estimate of the G-parameter: |G|/vr ~ 0.1 /a. This is about an order of magnitude 
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larger than is required for the charge separation [5]. Due to the electric field Eq of the parity-odd 
bubble quarkonium dissociation rate is finite at low P±, as indicated in Fig.[TJ (Fortunately, J/^'s 
can be measured down to very low P± [26\ [2"7] ) . The effect of the electric field is most clear along 
the direction perpendicular to the reaction plane, because magnetic component of the Lorentz force 
vanishes in this direction, see Fig. (3^a). 

Azimuthal distribution of dissociation rate is strongly asymmetric in external magnetic field. 
The second harmonic v% of the azimuthal distribution of survival probability is large and negative 
at low P±, while at high P± it is positive; zero of t>2 depends on the relation between the electric 
and magnetic fields. According to the preliminary experimental data J/^'s V2 is either small, about 
a few per cent, or zero |28U29| . Absence of such asymmetry in the experimental data may have two 
reasons, (i) Magnetic field is significantly weaker and short-lived than suggested in [HE], which is 
however at odds with the charge separation observations [8l-ll0|. (ii) Almost none of J/ip J s produced 
in the center of QGP survive. Rather they originate from the peripheral regions. The later scenario 
is realized if time t spent by quarkonium in the field is large because the dissociation rate only 
linearly increases with t but exponentially decreases with the decrease of the field strength toward 
the QGP periphery. Finally, if the bubble interference effects due to finite bubble size discussed at 
the end of Sec. [n] are important they can significantly reduce the CP-odd effect on J/ip dissociation 
rate. 

Abundance of possible effects associated with strong magnetic field calls for a detailed experi- 
mental investigation. 
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Appendix 



Here we list functions C and P that appear in nb\ 



12 



7_ 

TO 

exp 
f 1 



(to coth T + 

TO 



ln^ + 

. 2 7 



sinh to cosh to 

TO 

7 



- 2) sin 2 + sinh 2 t cos 2 9 
1 



-1/2 



T - T 



(t 2 -t 2 )^cos^ + t 2 



/ cosh to — cos t \ 2 / sinh t t 
\ sinh tq / \ sinh tq tq 



1/2 



sin 



(A.l) 
(A.2) 

(A.3) 



